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We consider random matrix models for the thermodynamic competition between chiral symmetry 
breaking and diquark condensation in QCD at finite temperature and finite baryon density. The 
models produce mean field phase diagrams whose topology depends solely on the global symmetries 
of the theory. We discuss the block structure of the interactions that is imposed by chiral, spin, and 
color degrees of freedom and comment on the treatment of density and temperature effects. Exten¬ 
sion of the coupling parameters to a larger class of theories allows us to investigate the robustness of 
the phase topology with respect to variations in the dynamics of the interactions. We briefly study 
the phase structure as a function of coupling parameters and the number of colors. 
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I. INTRODUCTION 

Some thirty years ago, it was observed that dense and cold quark matter might exhibit Cooper pairing as a result 
of an attractive quark-quark interaction in the color antitriplet channel [j]]. More recent models, based on non- 
perturbative effective interactions or on diagrammatic calculations of single-gluon exchange interactions, indicate 
that a color superconducting phase might develop pairing gaps as large as A ~ 100 MeV for quark chemical potentials 
on the order of 300 MeV Q|. This interesting possibility has direct consequences on the physics of dense stars and 
is certainly important for the determination of the phase diagram of nuclear matter under extreme (ultrarelativistic) 
conditions p, 0- 

Different order parameters have been proposed in the literature and studied as a function of, e.g., the quark masses, 
the number of flavors and colors, the quark chemical potentials, and temperature 11 In the limit of QCD with 
two flavors of light quarks (the 2SC limit), the order parameter has the form 

(^fcaip^gpcr 1 (~P)) = £ fg £ «/33 £ ctu' , (!) 

where p is a four-momentum and where we have displayed the flavor ( f,g ), color (a, /?), and spin (cr, a ) indices. 
The tensors £ ensures that the condensate is antisymmetric in flavor, color, and spin. Color is broken from SU( 3) to 
SU( 2), but the flavor symmetry SU(2)l x SU( 2)r and the baryon symmetry U(\)b remain intact. In the other limit 
of three degenerate flavors of light quarks, the favored order parameter exhibits a coupling between color and flavor 
rotations (called color-flavor locking or CFL) and has the approximate form §3 

('<Pfav(p)lp g p' r ' (~p)) = ®(p 2 ) £ fgA £ a 0A £ aa ' ■ (2) 

Color, flavor, and baryon symmetries are now broken down to SU(3) co \ oi +l+r x Z 2 . Since both SU(3)l and SU(3)r 
are locked to SU(3)color, the CFL condensate also breaks chiral symmetry. 

Following Ref. 0, the conjectured phase diagram for QCD with three flavors and realistic quark masses is given in 
Fig.m At asymptotically high densities, the scale is set by the quark chemical potential, p. For p m s ( m s is the 
strange quark mass), the more symmetric CFL phase is favored. As /i decreases, the increasing difference between u- 
(or d-) and s-quark Fermi momenta weakens (us) and (ds) condensates and eventually leads to a transition to the 
2SC phase @, possibly via a so-called LOFF phase characterized by a spatially varying gap i6:]. 

Most of these results have not been confirmed by lattice simulations, which are difficult to perform and interpret 
for QCD with three colors at finite baryon densities. The difficulty resides in the fact that the fermion determinant in 
the partition function is complex for non-zero p, and sampling weights are no longer positive definite. This difficult 
problem represents a significant barrier to understanding the phase structure of QCD in lattice simulations. 

Random matrix theory offers models that are capable of distinguishing those physical properties that are determined 
by global symmetries from those that depend on the detailed dynamics of the interactions. The random matrix 
Hamiltonian mimics the true interactions by adopting a block structure that is solely determined by the global 
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FIG. 1: Sketch of the phase diagram for three flavor QCD with realistic quark masses. 


After Ref. Qj. 


symmetries under consideration. The matrix elements are drawn on a random distribution, usually a Gaussian 
distribution, so that the model can be solved exactly. Such an approach leads to universal results that can be studied 
at two different levels. At the microscopic level, the statistical properties of the lowest eigenvalues of the Dirac 
operator are determined by the spontaneous breaking of chiral symmetry.as indicated for example by the universality 
of the spectral density near zero virtuality and related sum rules [7, St [13 • At the macroscopic level, many 
of the properties of the phase diagram (such as its topology and the presence of given critical lines or points) are 
independent of the detailed form of the interactions and are thus symmetry protected. The random matrix approach 
treats low-lying fermion excitations and, in its usual form, neglects their momentum and kinetic energy. The resulting 
phase dia gram is mean-field. Near critical regions, it produces results similar to those obtained in a Landau-Ginzburg 
approach [ 11L11 2l| . In general, random matrix models provide a useful tool for studying systems with non-trivial phase 
diagrams. 

This talk focuses on random matrix models that implement coexisting chiral and color symmetries in QCD with two 
light degenerate flavors (2SC). These models are extensions of the chiral random matrix models, which we introduce 
in Sec. im We present the construction of the color and spin block structure of the interactions in Sec. |111| The phase 
diagram is discussed in Sec. nvi This talk is a summary of previous work on the question, see Refs. [I3lll4l .[l5l|. 


II. CHIRAL RANDOM MATRIX MODELS 


We first consider chiral symmetry alone and turn to chiral random matrix models IBIS For QCD in the sector 
of zero topological charge with Nf flavors and zero chemical potential, the partition function is given as 


Z = 


N f 

DW Dip* Dip i exp 
2—1 


Nf 


2=1 




(3) 


where W is an n x n matrix which models the interactions. Its elements are drawn on a Gaussian distribution of 
mean zero and inverse variance E. Here, n is a measure of the number of low-lying degrees of freedom and is to be 
taken to infinity at the end of the calculations (i.e., in the thermodynamic limit). In the chiral limit, m = 0, the 
Dirac operator, D, has a block structure imposed by the chiral symmetry of QCD, { 11 , 75 } = 0. In the basis of the 
eigenstates of 75, this leads to the block structure 


V = 


0 iW \ 

iW* 0 )' 


(4) 


For random matrix models of QCD with SU (3) and fermions in the fundamental representation, W is complex. This 
choice corresponds to the chiral unitary ensemble, which is characterized by a Dyson index [3 = 2. The QCD Dirac 
operator in SU( 2) with fermions in the fundamental representation satisfies an additional anti-unitary symmetry, 


[C(a 2 ) color K 1 iV}=0, 


( 5 ) 
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with (C*(cj 2 )coior-^0^ = 1- (C is the charge conjugation operator, a 2 is the antisymmetric color matrix, and K is the 
complex conjugation operator.) This implies that it is possible to find a particular basis of states in which iV is 
real. Accordingly, W is chosen real in chiral random matrix models for SU( 2), and this leads to the chiral orthogonal 
ensemble with an index (3 = 1 E3- For fermions in the adjoint representation of the gauge group and any number of 
colors, the Dirac operator obeys the anti-unitary symmetry C~ l I\ with (C~ 1 K) 2 = —1. This leads to the symplectic 
chiral ensemble with (3 = 4 and quaternion real matrix elements MM- We will not consider this ensemble further 
in this talk. 

The essential difference among the ensembles lies in the number of independent random variables that are allowed 
per matrix element. (I.e., A complex number has two degrees of freedom; a real number has only one.) This difference 
is important in determining the statistical properties of the Dirac operator and, in turn, affects the phase diagram. 


III. RANDOM MATRIX MODELS WITH AN EXPANDED BLOCK STRUCTURE 

In order to study the competition between chiral, ( qq ), and diquark, (q T q) : condensates, we introduce an explicit 
dependence in the spin and color quantum numbers appearing in the 2SC order parameter of Eq. 0 We are thus 
lead to a Dirac operator with the chiral block structure of Eq. 0 , where W now has the following expanded color 
and spin sub-block structure: 


3 Nl-i 

W = EE A a (B> &n <8> A^ a . ( 6 ) 

fi—0 a= 1 

Here, the deterministic matrices represent spin and color degrees of freedom: a M = (1 ,ia) with a the Pauli matrices, 
whereas A Q are color matrices (Gell-Mann matrices for N c = 3). The random matrices, A^ a , are N x N and represent 
gluon fields. They are chosen real. Their elements are drawn on a Gaussian distribution with an inverse variance that 
is independent of /r and a in order to respect the Lorentz and SU(N c ) invariance in the vacuum. 

Including a quark chemical potential, /q a quark mass, to, and a temperature dependence, the partition function is 
now written as 


Z = f d 4 Dip *2 D0 2 t {[] DA^a} exp f -2IVS 2 ^ Tr[A, la A T ua ] 


x exp 


ipl 


iD + im + C+ 0 

0 —iD T — im — C'E 


if>i 

■02 


(7) 


where the dependence on the external parameters is given by C± = (=F<T 37 tT + ifx )74 with 


74 = 


0 1 
1 0 


( 8 ) 


Note that we have used a Gorkov representation and have transposed the flavor-2 fields so that the diquark condensate 
appears as an off-diagonal component of the inverse Dirac operator. 

Several comments are in order regarding the q and T dependences. The q dependence mimics the ni/j'-y 40 term 
in the Euclidean QCD Lagrangian [lj. The T dependence is based on two assumptions. First, temperature is 
introduced via a sum over the fermion Matsubara frequencies, uj n = imrT with n odd. As random matrix theory only 
treats low-lying degrees of freedom, much of the critical physics can be captured if one restricts the frequency sum 
to its two lowest terms 0 , hence the two-dimensional form a^irT = diag( 7 rT, —nT). The second assumption comes 
from the physical observation that the diquark order parameter in Eq. 0 couples fields of opposite four-momenta. 
We explicitly impose this coupling in Eq. 0 by taking opposite Matsubara frequencies for each flavor. Different 
temperature dependences have been proposed in closely related random matrix models m S3 ; we discuss their 
relationship to the present models in [ 2 l| ■ 


IV. PHASE DIAGRAMS 

The partition function can be evaluated following a classical method. An integration over A^ a produces a four- 
fermion interaction which can be Fierz-transformed to obtain the chiral and diquark channel terms. The resulting 
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FIG. 2: Phase diagram in the random matrix model for QCD with three colors and massless quarks. 


four-fermion potential can be expressed as a fermion field bilinear via a Hubbard-Stratonovitch transformation. These 
successive steps yield ® 

Z ~ j dadA e ~ ANn ^ A \ (9) 

where the thermodynamical potential, 0 (<t, A), is given as 

AT — 9 „_. __. 

fl(cr, A) = AA 2 + Bcr 2 -A-— ^ log ((a + m ± ^) 2 + tt 2 T 2 ) — ^ log ((a + m ± fi) 2 + n 2 T 2 + A 2 ) . (10) 

Here, A is the auxiliary field associated with the diquark channel (with condensates (V , 2 f{(* cr 2 )spin ^ 2 ^ir) = 

(V , 2 Z,(* cr 2 )s P in ^ 2 ^il)), while o' is related to the chiral channel (with condensates (VqV’i) = ^ 2 )) • The parti¬ 

tion function can be treated exactly in the thermodynamic limit N —> 00 by a saddle point method. The equilibrium 
values of the field then satisfy the two gap equations 


on 

= 0 , 

( 11 ) 

dVL 

dA 

= 0 , 

( 12 ) 


which constitute a system of polynomial equations that can be solved analytically or numerically. 

The form of the potential in Eq. m has a straightforward interpretation. The quadratic terms correspond to the 
energy cost for having constant auxiliary fields. The logarithmic terms are to be related to the single quasiparticle 
spectrum in a given set of fields er and A. Keeping in mind that the fermion four-momenta are neglected, these 
energies are given as e = ((<7 + m- F p ,) 2 + A 2 ) 1 / 2 for the two pairing colors (with the =p sign standing for either quark 
or antiquark excitations). For the N c — 2 unpaired colors, e = |cr + m =p fi\. 

We now turn to a chiral symmetric theory with m = 0. Note that, because it is possible to rescale all fields and 
external parameters by a constant, the topology of the phase diagram depends only on the ratio B/A. This is in 
fact a ratio of the Fierz constants in the chiral and diquark channel, and it measures the relative importance of the 
two symmetries. For example, small ratios are obtained with A^> B; the energy cost of AA 2 then prohibits diquark 
condensation. Similarly, large ratios B/A disfavor chiral symmetry breaking. For N c = 3, the interactions in Eq. ® 
lead to a ratio B/A = 3/4, which corresponds to the phase diagram of Fig. 0 One observes a low-density chiral 
broken symmetry phase separated from a higher density diquark phase by a first-order line. As the temperature is 
raised, the diquark phase makes a transition to the chiral symmetric phase across a second-order line. The transition 
from broken to restored chiral symmetry is second-order at low densities and first-order at intermediate densities. 
The two lines are separated by a tricritical point, in the vicinity of which the thermodynamical potential reduces to 
a (f> 6 theory m®. 

It is interesting to ask how the phase structure evolves if one changes the channel couplings A and B from the values 
representative of QCD. To this end, we have considered Hermitean Dirac I'D operators spanning an exhaustive set of 
combinations of helicity, spin, and color block structures. Remarkably, this set produces ratios B/A in the bounded 
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range [0,N c /2], This result differs from what would have been obtained in a pure Landau-Ginsburg approach, 
where there is no a priori knowledge of existing constraints among the coefficients of the effective thermodynamical 
potential. Here, however, having started at a more microscopic level, we are capable of discovering possible bounds 
on the coupling ratios. Each ratio corresponds to a separate phase structure. The major conclusions resulting from 
the study of phase diagrams for ratios in the allowed range are as follows: 

• there is only a finite number of different topologies. As B/A is varied continuously, the evolution from one 
topology to another is marked by the emergence or the vanishing of new critical points or lines; 

• it takes moderate — but finite — alterations of the theory to depart from the topology of Fig. El In that sense, 
the phase structure of Fig. Q is protected by symmetry. 


V. DISCUSSION 

A. Comparison with a microscopic theory 

In order to appreciate the approximations involved in the formulation of the random matrix model, consider the 
gap equation in the limit a = 0, 

30 2A 

<9A “ ^ “ A 2 + g 2 + 7r 2 T 2 ’ 

and compare with that obtained in a microscopic mean-field theory such as that of 
interaction modeled by that induced by instantons. Approximately, we find 

d A p f 2A 2A 

(2^)4 [a 2 + ( P ~p) 2 +pI + A 2 + (p + g) 2 +p 2 

where we have dropped the form factors needed for the convergence of the integral. It is also implied that the integral 
over p 4 is a sum over Matsubara frequencies, p^ = mrT with n odd. 

Because we have neglected the fermion four-momenta, the gap equation in the random matrix approach does not 
contain an integral. This has two consequences. First, the right term in Eq. m does not exhibit the logarithmic 
divergence observed in Eq. Q for p ~ p and A = 0. This divergence arises for values of p near the Fermi momentum 
in the limit A —> 0 and implies that a non-zero A must develop at all p. In contrast, the diquark phase in Fig. El 
does not exist for asymptotically high values of p. The random matrix interactions saturate in the diquark channel 
for p larger than ~ E. The second consequence of the absence of a momentum integral is that the gap, A, evolves 
monotonically as a function of p in contrast to results from a microscopic approach, see Refs. j2S| Ell ■ 

These discrepancies should not be considered as weaknesses of the random matrix approach. In the first case, A 
vanishes in a region where microscopic theories tend to produce small gaps which probably do not survive fluctuations 
beyond mean field. What random matrix cannot reproduce, however, is the A ~ e -1 /® behavior as a function of the 
QCD coupling constant, g. This behavior is due to the magnetic gluon interactions 1 2 ll . The random matrix approach 
is also unable to reproduce the unbounded increase of A at asymptotically high p , which is related to the running of 
9 EH- Both these features are predicted from diagrammatic theory and are related to particular dynamic processes. 
The second observed difference, the variation of gaps at moderate values of p, arises from the detailed dynamics of 
the interactions and depends on the choice of the regulating form factors. This behavior is not dictated by symmetry, 
and it is thus no surprise that it is not revealed by the random matrix approach. 


A ~ G 


(13) 

Ref. E3, based on an effective 

) , ( 14 ) 


B. QCD with two colors 

In the limit of N c = 2, the gauge group is pseudoreal, and quark and antiquark states transform similarly under 
global color rotations. They can be combined into spinors which obey an extended flavor symmetry SU (2Nf) for which 
(qq) mesons and (qq) baryons belong to the same multiplets. The random matrix model (with Nf = 2) reproduces 
this extended symmetry in the vacuum and its breaking pattern as a function of /i, T, and m. In the vacuum and in 
the chiral limit, the thermodynamic potential depends on the condensation fields through the combination a 2 + A 2 . 
The extended 517(4) symmetry is here apparent since a state with (a, A) = (E, 0) is indistinguishable from its rotated 
version with (cr, A) = (0,E). For m > 0 and low temperature, the phase (E,0) undergoes a second order phase 
transition to a diquark phase at p c — ra^ll where m„ ~ (mE) 1 / 2 is the pion mass |14| . Many results of the random 
matrix approach agree with chiral perturbation theory in this case l2(i . 
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C. Two more questions 

The random matrix model is a theory of low-lying modes. We argued earlier that restricting the sum over Matsubara 
frequencies leads to a model which nicely captures the critical physics. The neglect of high-energy modes leads however 
to unphysical results, such as a negative baryon density and a variation of the chiral field as a function /i, both in a 
theory with N c = 2 and in the region /i < m^/2. We have shown that the inclusion of appropriate high-energy terms 
(which should not describe the critical physics), either in the form of correction terms or as a sum over all Matsubara 
frequencies, fixes these anomalies while leaving the topology of the phase diagram intact 0 

Another question is whether the interactions in Eq. preserves the statistical properties of the eigenvalues of the 
Dirac operator that are expected from chiral symmetry alone. These properties should follow the predictions of the 
chiral unitary ensemble for QCD with three colors and those of the chiral orthogonal ensemble for QCD with two 
colors. This question is related to the number of random degrees of freedom allowed for the matrix elements. We 
have shown that, even if the interaction matrices are complex for both N c = 3 and N c = 2, the deterministic spin 
and color dependences lead to the properties that are expected for the spectrum of the Dirac operator 0 


VI. SUMMARY 

We have considered random matrix models of QCD with two flavors that are capable of developing chiral and diquark 
condensation. We have studied the symmetry breaking patterns as a function of temperature and quark chemical 
potential. The phase diagrams can be established exactly from an analytical evaluation of the partition function. 
The resulting thermodynamical potential has a straightforward interpretation in terms of elementary excitations. 
Upon arbitrary variations of the coupling constants in the two condensation channels, the phase structure evolves 
continuously but can only adopt a fixed set of topologies. 

This study shows that random matrix theory provides useful tools for studying systems with non-trivial phase 
diagrams and for distinguishing those properties that are protected by symmetry. 
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